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Abstract 

We present a continuum model for doped manganites which consist of two species of 
quantum spin 1/2 fermions interacting with classical spin fields. The phase structure 
at zero temperature turns out to be considerably rich: antiferromagnetic insulator, 
antiferromagnetic two band conducting, canted two band conducting, canted one band 
conducting and ferromagnetic one band conducting phases are identified, all of them 
being stable against phase separation. There are also regions in the phase diagram 
where phase separation occurs. 
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1 Introduction 



Doped manganites La\_ x A x MnOj, (A divalent) |IJ] are receiving quite a lot of both theo- 
retical 0, H, H, |, ||, [7], ||, H and experimental [|l(J attention lately. These materials show an 
interesting interplay between magnetism and conductivity with intrincated phase diagrams 
which are still controversial. 

In a cubic lattice the 3d orbitals of Mn split into a t2 9 triplet and an upper e g doublet. 
Due to the electronic repulsion and the Fermi statistics (Hund's rule) the three ti g levels are 
always single occupied forming a core S = 3/2 spin. The e g orbitals may be further splitted 



by a static Jahn- Teller distortion at small doping [III 



The above features are encoded in the so called double exchange models of different 
degrees of complexity. The simpler ones assume a strong Jahn- Teller distortion so that only 
the lower e g level is consider. Hence there is a single fermion field in each site, with a spin 
independent hopping term and a local interaction with the core spin [|, [5]. Core spins also 



interact among themselves with the usual Heisenberg term. Under certain assumptions |L2 
the interaction with the core spin can be traded for an angle dependent hopping term P, |I3|. 
The next level of complexity consist of taking into account the two e g levels |], |7|], and 
only very recently, the Jahn- Teller distortion has been incorporated dynamically by some 
authors M. 

It is the aim of this work to present a simple continuum model for doped manganites 
which also encodes the basic features above and, moreover, is exactly solvable for classical 
core spins. It produces a rich phase diagram which is in qualitative agreement with recent 
results and it shows, in addition, that stable canted phases exist. The main advantage with 
respect to previous approaches is that all the parameters of the material (lattice spacing, 
band curvature, Hund coupling, Heisenberg coupling and doping) combine into only two 
constants. This allows to present a two dimensional phase diagram which holds for a large 
amount of materials. 



2 The Model 

Cooperative phenomena are amenable of a field theoretical description. When the phe- 
nomena do not depend on the details of the microscopic system but only on its long wave 
length behaviour a continuum field theory description is appropriated. The field theoret- 
ical continuum model must contain the relevant degrees of freedom at long wavelengths, 
which depend on the particular systems and phenomena that are to be studied. In our case, 
these are doped manganites and their phase diagram at zero temperature. These systems 
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are known to undergo a number of phase transitions when the doping is increased. They 
are insulating antiferromagnets (AFI) at zero doping and become conducting ferromagnets 
{FC) at large enough doping. What happens between these two regimes is still controver- 
sial, though most of authors agree that the phase diagram is very rich and non-trivial. Early 
works on the subject suggested that an interesting intermediate conducting canted phase 
exists [[13] , but recent experimental |10[ and theoretical 0, f| results indicate that the canted 
phase appears to be unstable against phase separation. 

Theoretical work on the subject is based on variations of the double exchange models. 
The phase structure of the system is obtained from these models using certain simplifying 
assumptions (slave boson formalism ||, trial wave functions Q,. . . ) or extensive numerical 
simulations @, the scope of which is difficult to evaluate. We present below a continuum field 
theoretical model which, as we shall argue, contains the relevant long wavelength degrees 
of freedom of the system. Then our main assumption is going to be that the rich phase 
diagram of manganites can be understood from long wavelength physics only. As the model 
is exactly solvable, there are no further uncertainties due to uncontrolled approximations. 

Since we wish our model to include the well established AFI and FC phases, we need at 
least an AF order parameter field, a F order parameter field, and a, I — C order parameter 
field. For the AF and F order parameter fields we shall use Mi(i) and M2(x) the local 
magnetisations in the even and odd sublattices respectively. Both in the AF and F phases 
these local magnetisations are smoothly varying fields. In the AF phase M!(i)M 2 (i) ~ — 1 
whereas in the F phase Mi(i)M2(i) ~ 1. For the I — C order parameter one could think of 
introducing a single slowly varying spin 1/2 fermion field together with a chemical potential 
which regulates the doping. When the chemical potential is below the energy gap of the 
lowest spin state we have an I phase, when it overtakes this energy gap we have a one 
band C phase, and when it overtakes the energy gap of the highest spin state we have a 
two band C phase. However, a spin 1/2 field naturally couples to the local magnetisation, 
which changes abruptly from the even to the odd sublattice in the AF phase. Hence in this 
phase a single spin 1/2 field cannot be slowly varying over the system. We need at least 
two slowly varying spin 1/2 fermionic fields, ipi{x) which couples to the magnetisation in the 
even sublattice Mi(x) and fafe) which couples to the magnetisation in the odd sublattice 
M2(i). Since the conductivity is due to fermions moving from one sublattice to the other 
one a (spin independent) hopping term is introduced. The allowed values of the chemical 
potential will be limited by the physical condition that no conduction must exist when the 
hopping parameter vanishes. 

The model must be SU(2) spin invariant since the magnetic interactions emerge from the 
usual superexchange mechanism together with the Hund's rule. The space-time symmetries 



2 



of the underlying crystal must also be implemented and will be the only remain of the 
microscopic lattice structure. For simplicity we shall take a cubic lattice and comment later 
on the slight modifications that occur for other crystals. 

The lagrangian of the model reads 



C(x) 



ip\(x) 
+ 4(x) 



1 + ie)id + ^ + fi + J H ^rM 1 (x) 
2m 2 



d 2 <r 

1 + ie)ido + 7^ + V + Jh-M 2 {x) 



+ t (i>\(x)4> 2 (x) + i/>l(x)i/> 2 (x)) - J AF M 1 (x)M 2 (x). 



(2-1) 



The size of the parameters in the model are estimated by comparing them with the naive 
continuum limit of lattice double exchange models. For a cubic lattice we have 2m ~ l/a 2 t l , 
t ~ zt l , Jh ~ J\j and Jaf ~ zJ l AF /a 3 > 0, where a is the lattice spacing, z = 6 is 
the coordination number and the superscript I means the analogous lattice quantity. The 
fields ipi{x) may describe either electrons or holes. Since the conduction in actual doped 
manganites is due to holes, one should better figure out ipi(x) as hole annihilating fields. 
Recall that for holes Jh is negative whereas it is positive for electrons. This sign however is 
going to be irrelevant as far as the phase diagram is concerned. 

The lagrangian above is invariant under the following transformations: 



(i) Global SU (2) spin transformations, 

Mf{ x ) — B$M?(x), 



1,2) 



(2.2) 



(ii) Primitive translations, 



ipi(x) 
Mi (a) 



ip 2 (x) ip 2 (x) — >i>i(x) 

M 2 (x) M 2 (x) -^Mi(x), 



(2.3) 



(iii) Point group transformations, given by the group m3m 

A{x) — ► gM^x) 
M?(x) — RtiOM^x) 



(i = l,2), 



(2.4) 



when the point group transformation £ maps points in the same sublattice, and 



^i(x) — ► g^ 2 (i 1 x) i/j 2 (x) 
Mf(x) — RtiOMK^x) M%(x) 



RtiOM^x), 



(2.5) 
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when the transformation £ maps points of different sublattices. Anyway, the rotations 
g% and B%(£) can be absorbed by a SU{2) transformation and the change of sublattice 
in ( |2.5| ) by a primitive translation. Hence, in practice, we only have to care about the 
transformation of the coordinates. 



(iv) Time reversal, 



Mi(x) — > -Mi(Tx) 
where Tx = (— t, x). 



C = e~ ma 11 = -%o l , (2 = 1,2), (2.6) 



3 Effective Potential 

In order to find out how the ground state of the system changes as a function of the 
chemical potential, we shall calculate the effective potential and minimise it with respect 
to the order parameters Mi and M 2 . We shall assume that the ground state configuration 
corresponds to constant magnetisations both in the odd and even sublattices. Hence, the 
effective potential is to be minimised with respect to the angle 9 between Mi and M 2 only. 
We use y = cos(6 l /2). When y = 0, < y < 1 and y = 1 we have an antiferromagnetic, 
canted and ferromagnetic phase respectively. 

The effective potential is obtained by integrating out the fermion fields in the path 
integral, and it is formally given by 

V eff = JafM!M 2 + itr hgO/VT, (3.1) 

where 

^ _ ( (1 + ie)id + df/2m + p. + ^<tMi t 

\ t (1 + ie)id + df/2m + \i + f(rM 2 

and the trace is both on spin indices and space-time coordinates. VT is the volume of the 
space-time. 

If O has eigenvalues A n 

trlog6 = ^logA n . (3.3) 

n 

We have then to diagonalise the operator O. Since it contains only constant fields the 
diagonalisation with respect to the space-time is trivially attained by plane waves. The 
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diagonalisation with respect to the spin indices is a simple linear algebra problem. We 
obtain 

k 2 

An = Oi(q) = {l + ie)u- — -Qi (3.4) 

2m 



«, = ± ^^ + ^ ± 2 7 cos|-, , 7 .^. (3.5) 

g = (o>,k) and M = |Mi| = IM2I = 3/2. The restriction for the values of the chemical 
potential in the model implies that at most the two lower eigenvalues in ( |3.5|) may contribute. 
This motivates the following reparametrisation of the chemical potential: 



\Jh\M 



y/l + 7 2 - 2 1Vo (-Ky <yr x = l/2), (3.6) 

which eases comparison with the energy levels in (|3.5| ) (y = cos(#/2)). In order to simplify 
the analysis we assume 7 small and keep only linear terms in 7 in the relevant eigenvalues 
above. Namely, 

a = -i^j^7(lfo±v). (3-7) 
This is justified for t <C Jh, as it turns out to be the case for the actual materials [fL4}] . Any- 



way, this simplification can be lifted with the only drawback that the few analytic expressions 
below must also be substituted by numerical analysis. 

In order to calculate the sum ( |3.3| ) we have used (^-function techniques which are 
explained in the appendix. We obtain the effective potential (for // < 0) 



V,n = Vr 



(2y 2 - 1) - A ((y + yf' 2 d{y + y) + (y - y f 2 9(y - y))] , (3.8 



where we have defined 



(2m) 3/2 t 5/2 z 3/2 t 

V « = Ja ? M2 ' A = 15n 2 J ApM 2= { J AFa 3 M iy ^ 



4 Phase Structure 



The possible phases of the model are obtained by minimising (|3.8| ) with respect to y 
for the different values of the parameters A and y^. The number of conducting bands is 
given by the number of ^-functions in ( |3.8j ) which contribute to the effective potential at the 
minimum. 

In order to gain some qualitative understanding and to make the minimisation procedure 
systematic we shall first separate the cases y < and y > 0. For each case we shall work 
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out the stability conditions for AF (y = 0), canted (0 < y < 1) and F (x — 1) phases. After 
that we shall compare the energy of the stable phases and obtain the curves which separate 
them. 

The stability conditions are given for the different phases by 

AF: V^(0)>0 or V^(0) = V% f (0) > 
C : V^(y c ) = K//(?/c) > (4.1) 
F : 1) < 0. 

Let us then consider first the case yo < 0. Clearly for y < —1 the unique existing phase 
is the AFI phase. In the case — 1 < yo < only the lowest of the four spin eigenvalues 
may contribute to the effective potential. The stability conditions yield the following stable 
phases: 

AFI : y = 

FC: y = l A(l + y f 2 >8/5. (4.2) 
The canted phase is not stable as it can be seen from the condition V^(y c ) = 0, 

y c = h -A(y Q + y c f'\ (4.3) 

which has at most one solution y c G [—yo, 1]. Since V e g is continuous, and increasing at 
y = this solution must be a maximum when it exists. 

The curve V e ff(0) = V e jj(l) in the plain (y ,A), which separates the AF and F phases, 
reads 

A(l + y,f' 2 = 2 (-Ky o <0). (4.4) 
Above this curve the F phase is favoured against the AF phase and viceversa. 
Consider next the case < yo < 1. The stability conditions are given by 

AFC2 : y = Ay X J 2 < 8/15 

CC2: 5A(y 2 + 3yl)/4=(y + y c r/ 2 + ( yo -y c )^ 8/15 < Ay]/ 2 < 2^/5 

CCl: y c = 5A(y + y c ) 3 / 2 /8 Ay/ 2 > 2^/5 (4.5) 

FC1: y = l A(l+y ) 3 / 2 >8/5, 

where AFC2, CC2, CCl and FCl stand for antiferromagnetic two band conducting, canted 
two band conducting, canted one band conducting and ferromagnetic one band conducting 
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respectively. Notice that AF and canted phases do not compete among them, but only with 
the F phase. The curves providing the boundary between the different phases are given by 

A[(l + y f/ 2 - 2y 5 /2 ] = 2 



AFC2 - FCl : A[{1 + y ) b/2 - 2y^1 =2 < y < 0.127195 

AFC2 - CC2 : Ayl' 2 = 8/15 0.127195 < y < 1 

CC2 - FCl : 5A(y 2 + 3y 2 )/4 = (y + y 2 f' 2 + (y - y 2 fl 2 0.127195 < y < 0.168457 

CC2 - CCl : Ayl 12 = 2^2/5 0.168457 < y < 1 (4.6) 

CCl - FCl : 5A(y + ?/i) 3/2 /8 = Vi 0.168457 < y < 0.5 

CCl - FCl : 5A{1 + y ) 3/2 /8 = 1 0.5 < y < 1, 

where y± and ?/2 are given implicitly by the equations 

[(i + y ) 5/2 - (yo + V2? 12 - (y - y 2 ) 5/2 ][(2/o + y 2 ) 3/2 + (y - y 2 ) 3/2 ] = |(i - y 2 2 )(y 2 2 + 3y 2 ) 

(4.7) 

(yx + 2/o) 5/2 + 2(1 + y ) 1/2 (yi + yo? + 3(1 - y Q )( yi + y ) 3/2 

+ 4(1 - 2y Q )(l + y ) 1/2 (yi + y ) - 8y (l + y )(yi + yo) 1/2 - 4y (l + Z/o) 3/2 = 0. 

The outcome is plotted in fig. [I]. 

Recall that fig. [I] actually does not plot a phase diagram against doping but against yo 
which is related to the chemical potential rather than to the number of conducting fermions 
or doping. Recall also that V e g is to be regarded as a (zero temperature) grand canonical 
potential rather than as a free energy. The doping is introduced via 



x = —a 



djj, 



a 3 dV ( 



<-ff 



t dy 



provided that one molecule exists per unit cell with a lattice parameter a. Taking into 
account ( |3.9|) the doping corresponding to the different phases reads 



AFI : 


x = 









z 3/2 


AFC2 : 


x = 






6tt 2 






z 3/2 


CC2 : 


X = 






6tt 2 






z 3/2 


CCl : 


X = 






6n 2 






z 3/2 


FCl : 


X = 






6tt 2 



iyt 



[(yo + y c ) 3/2 + (yo - y c ) 3/2 } 
(yo + y c f 12 
(i + yof 12 . 



(4.9) 



where the y c for the CC2 and CCl phases are given in ( (4.5[ ). 
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Figure 1. Phase diagram in the (yo,A) plane. The thick solid line corresponds to first order 
transitions whereas the remaining solid lines to second order ones. The dotted and dashed dot- 
ted lines are the upper stability boundaries for the CC1 and CC2 phases respectively. The two 
dashed lines are the boundaries for the reliability of our model for z\ Jh\M/2( Jafo? M 2 ) ~ 50 and 
z\ Jh\M/2( Jafo, 3 M 2 ) 200 respectively. Only the part of the phase diagram to the left of the 
corresponding dashed line is trustworthy in each case. 

These expressions for the doping permit us to establish that all our phases are thermo- 
dynamically stable, unlike the ones observed in ref. [§, [§]. This is easily proven from the 
stability condition dfi/dx > 0. For the F and AF phases this is trivially obtained, whereas 
canted phases are stable if they are below the curves 

CC2: 5Ay/3=(y + yf 2 -(y -y) 1 / 2 

y 2 - 5y 2 + 4y (y 2 - y 2 ) 1 ' 2 = (y < y ) (4.10) 
GC\ : AyJ 2 = 16/15V3. 

This is always the it is shown in fig. [I] where we have plotted the two curves. 

Once we have the expressions Q4.9Q for the doping it is straightforward to translate fig. [I] 
to a more conventional phase diagram where the doping, x, appears in one of the axes. This 
is given in fig. |^ (recall z = 6). 



Figure 2. Phase diagram in the (x, A) plane. PSi (i = 1, 2, 3, 4) indicates the new regions where 
the phases at their boundary may coexist. The x = axis corresponds to the AFI phase. The 
two dashed lines are the boundaries for the reliability of our model for z\ Jh\M/2( Jafo^M 2 ) ~ 50 
and z\Jh\M /2(Jafcl 3 M 2 ) ~ 200 respectively. Only the part of the phase diagram to the left of the 
corresponding dashed line is trustworthy in each case. 

It is interesting to notice that in fig. |] new regions arise, which we have denoted PSi 
(i = 1,2,3,4), between the FC1 and the AFI, AFC2, CC2 and CC1 phases respectively. 
This is due to the fact that the thick solid line separating FC1 and AFI, AFC2, CC2 
and CC1 in fig. |l| corresponds to a first order phase transition. Along this line two stable 
inequivalent minima have the same energy and the chemical potential cannot be traded by 
the doping. These regions are likely to consist of coexisting domains where the two phases 
at the boundary are realised (phase separation) IJ. AFI and FC1 would coexist in PSI, as 
it has been observed in recent works || £§. FC1 and AFC2, CC2 and CC1 would coexist 
in PSI, PS2 and PS3 respectively. These three last possibilities of phase separation have 
not been found before. 

As mentioned in section ^|, the fact that for t = we do not permit conductivity restricts 
the values that the chemical potential takes to y < y™ ax = 7/2. By substituting this 
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expression in A we obtain 

A - ^ zIJhIM (A 11) 

A ~ 15tt* 2(J AF a 3 M 2 ) Vo ■ (4 - U) 

which gives the boundary of validity for our results. It turns out to be a straight line in fig. [1] 

provided that Jaf and Jh remains constant as y™ ax moves, which can be straightforwardly 

translated to fig. |2|. Only the phase diagram to the left of this curve is trustworthy. 

We take for the coupling constants t/(J AF a 3 M 2 ) ~ 10 - 20 and z\J H \M /2(J AF a 3 M 2 ) ~ 
50 — 200, which is compatible with the values given in the literature. For these values 
A ~ 1 — 2, and the two extreme validity curves are displayed as dashed lines in fig. [I] and 
fig- 0. 



5 Conclusions 

We have presented a simple model in the continuum which is able to describe the rich 
phase structure of doped manganites for a wide range of these materials. 

We have assumed an underlying cubic crystal for simplicity. Nevertheless, the ortho- 
rhombic distortion can be easily accommodated by the following simple changes in the phys- 
ical parameters: m 3 — > m x m y m z , a 3 — > abc, J AF Jx + Jy + Jz and t — > t x + t y + t z . In 
practice this does not modify our results since it would only lead to a different A, which is 
anyway a free parameter in our phase diagrams. This fact also suggests that the structural 
transitions that these materials undergo when increasing the doping [FJ] are not essential in 
order to understand the F — AF and I — C transitions. 

An important feature of our results is that the two canted phases that we observe are 
stable against phase separation, unlike in some previous works |J. We also observe 
regions in the phase diagram where phase separations of several kinds may occur. If we 
plug realistic values for the physical parameters we find A ~ 1 — 2. Within this range the 
following sequences of phases are possible upon increasing x: (i) AFI — PS1 — FC1, (ii) 
AFI - AFC1 - PS2 - FC1, (iii) AFI - AFC2 - CC2 - PS3 - FC1, (iv) AFI - AFC2 - 
CC2 — CC1 — PSA — FC1. Recall also that in PS3 and PSA ferromagnetic and canted 
phases coexist, which is a situation that has not been contemplated in previous works. This 
may explain some controversial results obtained by different authors. 

Let us also mention that the two fermion fields ipi(x) and ip2{%) accommodate the e g 
doublet in our model. Indeed in the AF phase the two lower and two higher eigenvalues 
( |OD are degenerated. In the F and C phases the degeneracy is lifted. This implies that 
the splitting between the two e g levels receives a contribution from the dynamics of the 
conducting fermions in addition to that from the static Jahn- Teller distortion. 
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The model can be used in the future to study the temperature dependence of the phase 
diagram. Fluctuations due to spin waves in all the phases (including the canted ones) can 
also be incorporated |fT6| . It would also be interesting to see if the model can be generalised 
to accommodate the Jahn- Teller distortion dynamically. 
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Appendix A: ^-function Techniques 

The ^-function techniques provide a very efficient way to calculate the trace of the loga- 



rithm of operators [15|. The ^-function associated to an operator O is defined as 



(a.i; 



Then 



d 



is 



(A.2) 



Consider the operator O in ( |3.2| ). Once the spin diagonalisation is performed we only 
have to consider the space-time trace over a generic spin eigenvalue denoted by Oj. Since 
the real part of the operator — iOi is positive for positive energies and negative for negative 
ones, due to the term ieu in (|3.4|), it is convenient to consider the integral form of (q(s) over 
positive and negative energies separately. 



*r[(0*0(-"))~ 

tr[{6ie(u)y 



-t 



, „ , f° dw d 3 k 
drr / — — 

T(s) Jo J~oc 2vr (2tt) 3 



r s) Jo 



00 , 8 _i f°° dw d 3 k 



2ir (2vr) 3 



e -iOM)r y T ( A-3a ) 



After the energy and momentum integration we obtain the expressions 



tr[(diO(-u))- 

tr[{6ie(u)y 



VT /2m\ 3 / 2 r(s-5/2) 



167T \ 71 ) 



T(s) 



VT /2m\ 3 / 2 T(s-5/2) 



16vr V 



7T / 



r s) 



(-^) s+5 / 2 (^ l )- s+5 / 2 . 



(A.4a) 
(A.4b) 
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We need the derivative of the above with respect to s at s = 0. The presence of 1/T(s) ~ s 
makes the evaluation very easy, giving rise to 



d 
ds 



- - CdW 



s=0 



VT(2m) 
30^ 



3/2 



i5/2(_ ini )5/2 _ (_i)S/2( in .)5/2- 



(A.5) 



The expression between square brackets vanishes when f2, > 0, i.e., when the chemical 
potential is bellow the energy of the z-th state, and is non zero when Qi < 0, i.e., when the 
chemical potential is above the energy of the z-th state. This leads to the effective potential 
(for n<0,y o < y™ ax < 1) 



V P , 



Va 



A 



7 



5/2 



'1 + 



_2jy_ 

1 + 7 2 



271/q 
1 + 7 2 



5/2 



0(yo + y) 



272/ 



"4 1 



272/0 



5/2 



9{yo-y) 



(A.6) 



VV 1 + 7 2 V ! + 7 2 
where y = cos(8/2), whereas 7, yo, and Vb are defined in ( ^5|) , flOl) and ( p.9|) respectively. 

(2m) 3 / 2 



.4 



15vr 2 J AF M 2 



5/2 



(A.7) 



Eq. ( p.7|) follows from the above by keeping only terms linear in 7. 
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